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Abstract. I summarize two recent applications of chiral perturbation theory to 
electromagnetic reactions on deuterium: elastic electron-deuteron scattering, and 
Compton scattering on deuterium. Both calculations have now been carried out to 
three orders in the chiral expansion. The expansion shows good convergence and is 
able to reproduce data for |q| < 600 MeV in ed and for u> = 55- 95 MeV in 'yd. These 
results demonstrate that yPT can be used to reliably compute operators and wave 
functions for low-momentum-transfer reactions in light nuclear systems. 


1. Introduction 

In recent years effective field theory (EFT) techniques have found increasing acceptance 
in nuclear physics. EFT facilitates the systematic separation of the effects of high- 
momentum physics from those of low-momentum physics. In nuclear physics there are 
many different momentum scales present. But for the A = 2 system the key ones would 
seem to be: 

VMB <C rrin A xSB , (1) 

with B the binding energy of deuterium, m n and M the pion and nucleon masses, 
and A x sb the scale of chiral-symmetry breaking, which is of the order of the mass 
of the p meson. An EFT based on the first hierarchy—V MB -C m n —has had 
considerable success describing very-low-energy reactions in A = 2 and A = 3 (for 
reviews see Dim and has recently been extended to A = 4 |3j. Here I will focus on 
the EFT built on the second hierarchy— m n -C A x g B . This EFT is chiral perturbation 
theory (%PT) (for a review see |S|). Here the low-momentum physics is described by a 
quantum field theory of nucleons and pions whose Lagrangian is the most general one 
that respects the approximate chiral symmetry of QCD and the pattern of its breaking. 
Higher-energy effects of QCD appear in yPT as non-renormalizable contact operators. 
Any process involving nucleons and an arbitrary number of (soft) pions and photons 
should then be computable as an expansion in the ratio of nucleon or probe momenta 
(denoted here by p and q) and to A x g B . 

Since typical momenta in light nuclei are < we would expect to be able to 
calculate the response of such nuclei to low-momentum electroweak probes using yPT. 
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As we will see below, %PT provides a systematic way to compute what in nuclear- 
structure physics would be called “bare operators”. But the operators computed in 
this fashion include pieces which account for the fact that they are not actually “bare” 
if one thinks about the physics at scale A x sb which has been “integrated out” of the 
low-momentum theory of nucleons and pions: high-momentum nucleon modes, the delta 
isobar, mesons such as the rho, pomerons, M-branes, etc. 

Here I will argue that these operators produce results for electromagnetic few- 
nucleon reactions which are model independent, can be systematically improved by 
going to higher order, and yield a good description of the data over the region of validity 
of the EFT. Section 2 outlines the yPT expansion, and briefly discusses how to obtain 
deuteron wave functions using it. Section 3 then looks at electron-deuteron scattering 
in yPT as a probe of deuteron structure. Section 4 examines the use of Compton 
scattering on the deuteron—in particular as a way to extract isoscalar polarizabilities 
of the nucleon. 

2. Power counting and deuteron wave functions 

2.1. Power counting 

Consider an elastic scattering process on the deuteron. If O is the transition operator 
for this reaction then the amplitude in question is simply ('if\0\if) 1 with \i[}) the deuteron 
wave function. In this section, we follow Weinberg J|, and divide the formulation of a 
systematic expansion for this amplitude into two parts: the expansion for O, and the 
construction of \if). 

Chiral perturbation theory gives a systematic expansion for O of the form 

OO 

0 = EO ( ” ) . (2) 

71—0 

where we have labeled the contributions to O by their order n in the small parameter 
P = p/ A x sb, <?/A x sb, "Itt/A x sb- 

To construct (PP one first writes down the vertices appearing in the chiral 
Lagrangian up to order n. One then draws all of the two-body, two-nucleon-irreducible, 
Feynman graphs for the process of interest which are of chiral order P n . The rules for 
calculating the chiral order of a particular graph are: 

• Each nucleon propagator scales like 1/P (provided that the energy flowing through 
the nucleon line is ~ m f ); 

• Each loop contributes P 4 ; 

• Graphs in which both particles participate in the reaction acquire a factor of P 3 ; 

• Each pion propagator scales like 1/P 2 ; 

• Each vertex from the nth-order piece of the chiral Lagrangian contributes P". 

In consequence more complicated contributions involving two-body mechanisms, and/or 
higher-order vertices, and/or more loops are suppressed by powers of P. 
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This leaves us with the problem of constructing a deuteron wave function that is 
consistent with the operator O. Weinberg’s proposal was to use the yPT expansion m 
for the NN potential V, and then solve the Schrodinger equation to find the deuteron 
(or other nuclear) wave function p]. Matrix elements fip\0\i]f) should then incorporate 
the consequences of chiral-symmetry breaking in a model-independent way. 

Calculations have shown that the NN phase shifts can be understood, and deuteron 
bound-state static properties reliably computed, with NN potentials derived from 
yPT |7, IHJ 0 M IB- This is now a sophisticated enterprise, with N 3 LO [0(P 4 )] 
potentials recently having been obtained du n2i nni- These potentials can reproduce 
NN data with an accuracy which rivals that of “high-quality” NN potentials, at least 
for Ti a b < 200 MeV. 

In the Goldstone-boson and single-nucleon sector yPT works because loop effects 
are generically suppressed by powers of the small parameter P. In zero and one-nucleon 
reactions the power counting in P applies to the amplitude, and not to the two-particle 
potential. However, the existence of nuclei tells us immediately that a power counting 
in which all loop effects are suppressed cannot be correct for the multi-nucleon case, 
since if it were there would be no nuclei, and neither we, or this workshop, would 
occur. Weinberg’s proposal to instead power-count the potential is one response to 
this dilemma. However, its consistency has been vigorously debated in the literature 
(see (DM for reviews). Beane et al. [II] have shown that Weinberg’s proposal is 
consistent to leading order in the 3 S i — 3 D\ channel, but there remain questions about 
the consistency in other channels and at higher orders ,0 US]. Here we will use wave 
functions calculated according to the Weinberg proposal. One of our objectives will be 
to assess whether such a procedure really does lead to convergent EFT expansions for 
the processes under consideration, or whether problems along the lines of those raised 
in Refs. rm mn also occur in the chiral EFT for electroweak reactions on deuterium. 


3. Elastic electron scattering on deuterium 


One quantitative test of yPT’s deuteron wave functions is provided by clastic electron- 
deuteron scattering. We thus turn our attention to the deuteron electromagnetic form 
factors Gc, Gq, and Gm 


These are matrix elements of the deuteron current J^, with: 


Gc = 

r< _ 1 

^Q 2\e\r)M'j 


b(0 

«o 


/ 


jjl) + (o 
°>-< 1 




J c 


o) + ( 

0). 


-1 




-1 


Gm — 


V^M 
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( 3 ) 

( 4 ) 


where we have labeled these (non-rclativistic) deuteron states by the projection of the 
deuteron spin along the direction of the momentum transfer q and r] = |q| 2 /(4Mj). 
Gc, Gq, and Gm are related to the experimentally-measured A, B, and T 2 o in the usual 
way, with T 20 being primarily sensitive to Gq/Gc and B depending only to Gm- Here 
we will compare calculations of the charge and quadrupole form factor with the recent 
extractions of Gc and Gq from data jT7j. 
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Both of these form factors involve the zeroth-component of the deuteron four- 
current J°. Here we split J° into two pieces: a one-body part, and a two-body 
part. The one-body part of J° begins at order \e\ (the proton charge) with the 
impulse approximation diagram calculated with the non-relativistic single-nucleon 
charge operator for structureless nucleons. Corrections to the single-nucleon charge 
operator from relativistic effects and nucleon sub-structure are suppressed by two powers 
of P, and thus arise at 0(eP 2 ), which is the next-to-leading order (NLO) for Gc and 
Gq. At this order one might also expect meson-exchange current (MEC) contributions. 
However, all MECs constructed with vertices from C^ N are isovector. These play a role 
in, e.g. np —> dy na- The first MEC effect in ed —> ed does not occur until N 2 LO, 
or 0(eP 3 ), where an NNn'y vertex from Pjfy gets included in an isoscalar two-body 
contribution to the deuteron charge operator (see Fig. mb)) t. 



Figure 1 . Diagrams representing the leading contribution to Jo [(a)], the leading 
two-body contribution [(b)], and the dominant short-distance piece [(c)]. Solid circles 
are vertices from and the shaded circle is the vertex from df'f. 


The most important correction that arises at NLO is the inclusion of nucleon sub¬ 
structure. At 0(eP 2 ) the isoscalar nucleon form factors are dominated by short-distance 
physics, and so the only correction to the point-like leading-order result comes from 
the inclusion of the nucleon’s electric radius. For the isoscalar combination of nucleon 
electric form factors \PT to 0(eP 2 ) gives: 

G e \PT NLO = 1 “ |( r E )2 )^ 2 - ( 5 ) 

This description of nucleon structure breaks down at momentum transfers q of order 

300 MeV. There is a concomitant failure in the description of ed scattering data mm- 

In order to focus on deuteron structure, in the results presented below I have chosen to 

circumvent this issue by using a “factorized” inclusion of nucleon structure t'SJ: %PT is 

used to compute the ratio -%j-. This allows us to use experimentally-measured single- 

Pe 

nucleon form factors § in the calculation, thereby allowing us to test how far the theory 
is able to describe the NN dynamics. 

f This exchange-charge contribution was first derived by Riska [111) . 

§ There is a bit of an issue of circularity here, since ed scattering data is one input to the extraction 
of the neutron electric form factor. 
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-LO: 0(e) 



Q (MeV) 


Figure 2. The deuteron charge and quadrupole form factors to NNLO in chiral 
perturbation theory. The experimental data is taken from the compilation of Ref. PH 
Gq is in units of fm 2 . 


The results for Gc and Gq are shown in Fig. [21 The figure demonstrates that 
convergence is quite good below q ~ 700 MeV. The results shown are for the NLO 
chiral wave function. It is clear that—provided information from single-nucleon form 
factors is taken into account—yPT is perfectly capable of describing the charge and 
quadrupole form factors of deuterium at least as far as the minimum in Gq- 

Gm can be obtained in a similar way. The LO contribution to Gm is O(eP). 
The first two-body mechanisms enter at 0(eP 4 ), when an undetermined two-body 
counterterm, as well as a long-range two-body current with an undetermined coefficient, 
appear m 12a. Results for Gm at 0(eP i ) turn out to be of similar quality to those for 
Gc EU, but are somewhat more sensitive to short-distance physics, as expected given 
the presence of a counterterm at NNLO in this observable. 

Using the NNLO chiral wave function, high-quality potential-model wave functions, 
or indeed simple wave functions which include only one-pion exchange regulated by a 
square well at short distances, does not modify the plots for Gc and Gq greatly below 
q = 700 MeV [21]. As long as wave functions with the same long-distance parts are 
employed these observables are apparently not very sensitive to the choice of | ip). In the 
language of the no-core shell model they do not “renormalize strongly”. (An important 
exception is the deuteron quadrupole moment, Q d .) This is good news for Weinberg’s 
power counting, because the first operators with which Jq can be renormalized occur 
at 0(eP 5 ) (Fig. 12(c)). Therefore it was to be expected that operator renormalization 
would have this small impact on Gq and Gq. Furthermore, the dominant piece of the 
operator renormalization can be taken into account by adjusting the coefficients of the 
0(eP 5 ) operators to reproduce, e.g. (r d ) and Q d - This would then push the omitted 
effects of operator renormalization in these observables to 0(eP 7 ). 
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Compton scattering on the nucleon at low energies is a fundamental probe of the long¬ 
distance structure of these hadrons. This process has been studied in \PT in 
The following results for the proton polarizabilities were obtained at 0(e 2 P ): 


CM.p 


5 e 2 g\ 


= 12.2 x l(r 4 fnt 


B v = — a v = 1.2 x 10 
xp 10 p 


-4 


fm 4 


( 6 ) 


3847 r 2 ffm 

Going to next order in the chiral expansion single-nucleon counterterms which shift the 
polarizabilities enter the calculation. Using the y p amplitude computed by McGovern 
to 0(e 2 P 2 ) [25] we fitted yp data over the kinematic range ui, \J\t\ < 180 MeV 


d(T/dll | ab 



^lab 

Figure 3. Results of the 0(e 2 P 2 ) EFT best fit to the differential cross sections for 
Compton scattering on the proton at various angles, compared to the experimental 
data from various facilities. The gray region is excluded from the fit. References for 
the data and a legend for the symbols are given in Ref. ESt- 

The best fit for the proton electric and magnetic polarizabilities is 

= (12.1 ± 1.1)1°;® x 10~ 4 fm 3 , f3 p = (3.4 ± 1.1)1°;* x 10 -4 fm 3 . (7) 

Statistical (1-er) errors are inside the brackets. An estimate of the theoretical error due 
to truncation of the expansion at 0(e 2 P 2 ) is given outside the brackets. The result (0) 
is compatible with other extractions, although the central value of /3 p is higher EH- 
One would like to perform a similar analysis for yd scattering data. Coherent 
Compton scattering on a deuteron target has been measured at E 1 = 49 and 69 MeV 
by the Illinois group |2B], at E 1 = 84.2 — 104.5 MeV at Saskatoon [22], and at E 1 = 55 
and 66 MeV at MAXlab at Lund ra- Such data should permit an extraction of 
the isoscalar nucleon polarizabilities ajv and /hv- These are interesting, not only as 
fundamental nucleon-structure parameters, but also because yPT predicts ajy = a p , 
Pn — (3 P at 0(e 2 P). The amplitude for yd scattering involves mechanisms other than 
Compton scattering on the individual constituent nucleons. So, our desire to extract 
nucleon polarizabilities argues for a theoretical calculation of this reaction that is under 
control in the sense that it accounts for all mechanisms to a given order in \PT. 
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The Compton amplitude we wish to evaluate is (in the 7 d center-of-mass frame): 

«(hh= T^(k',k) 

+ /if^ <M?') ryy, x (C£) v-m(p) (8) 

where M (.M') is the initial (final) deuteron spin state, and A (A') is the initial (final) 
photon polarization state, and k {k') the initial (final) photon three-momentum, which 
are constrained to\k\ = \k'\ = uj. The first integral is represented by the graphs of Fig.0 
where the photon interacts with only one nucleon. The second integral corresponds to 
the graphs of Fig. 0 where the current is of a two-body nature. 




I H H I 




Figure 4. Graphs which contribute to Compton scattering on the deuteron at 0{e 2 ) 
(a) and 0(e 2 P) (b-e). The sliced blob in graph (c) is from dff. Crossed graphs are 
not shown. 


The LO contribution to Compton scattering on the deuteron is shown in Fig. 0(a). 
This graph involves a vertex from and so is 0(e 2 ). This contribution is simply the 
Thomson term for scattering on the proton. There is thus no sensitivity to either two- 
body contributions or nucleon polarizabilities at this order. At 0(e 2 P) there are several 
more graphs with a spectator nucleon (Figs. 0(b), (c),(d)), as well as graphs involving an 
exchanged pion with leading order vertices (Fig. [3J) and one-loop graphs with a spectator 
nucleon (Fig. 0(e)) |TTT] . Graphs such as Fig. 0(e) contain the physics of the proton and 
neutron polarizabilities at 0(e 2 P) in yPT. 

We employed a variety of wave functions \ip), and found moderate wave-function 
sensitivity. Results shown here are generated with the NLO chiral wave function of 
Ref. [2J. Fig. 0 shows the results of a fit to the world’s modern yd data. The dotted line 
is the prediction at 0(e 2 ) in the kernel, where the second contribution in Eq. (JBJ) is zero, 
and the single-scattering contribution is given solely by Fig. 0(a). It is quite remarkable 
how well the 0(e 2 ) calculation reproduces the 49 MeV data. However, the agreement 
is somewhat serendipitous: there are significant 0(e 2 P) corrections. At these lower 
photon energies Weinberg power counting begins to break down, since it is designed for 
oj ~ tUtt, and does not recover the deuteron Thomson amplitude as u —> 0. Correcting 
the power counting to remedy this difficulty appears to improve the description of the 
49 MeV data, without significantly modifying the higher-energy results (2HIEU- 
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Figure 5. Two-body graphs which contribute to yd scattering at 0(e 2 P). Crossed 
graphs are not shown. 



Figure 6. Results of the 0(e 2 ) (dotted), 0(e 2 P) (dashed), and 0(e 2 P 2 ) (solid) 
calculations for yd scattering as compared to data at E 7 = 49 MeV J2S] ■ 55 MeV ;3D;, 
66 & 69 MeV | 28lld()| . and 95 MeV |29| . Only statistical errors on the data are shown. 
The gray region is excluded from the fit. 


At 0(e 2 P) all contributions to the kernel are fixed in terms of known pion and 
nucleon parameters, so to this order xPT makes predictions for deuteron Compton 
scattering. Note that the 0(e 2 P) corrections get larger as u is increased—as expected. 
The agreement of the 0(e 2 P) calculation with the intermediate-energy data sets is 
very good, although only limited conclusions can be drawn. These results are not very 
different from other, potential-model, calculations E2 ESI El- 

At 0(e 2 P 2 ) terms enter the single-nucleon amplitude which allow us to fit a n and 
/hv to yd data. A number of new two-body currents also appear. However, there are still 
no short-distance jNN operators contributing to yd scattering at this order. Therefore 
an 0(e 2 P 2 ) calculation allows us to test the single-nucleon physics which is used to 
predict the results of coherent scattering on deuterium, since there are no undetermined 
parameters in the two-body mechanisms that enter. 
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Unfortunately the data is fairly sparse, and the sensitivity to the NN wave function 
quite large, so our fit does not tie down a n and /3n very accurately. The result is: 

a N = (13.0 ± 1.9)1?;!! x 10 ” 4 fm 3 , Pn = (- 1.8 ± 1.9)1^ x 10“ 4 fm 3 . (9) 

The errors inside the brackets are statistical and those outside are an estimate of the 
effect from higher-order terms. (For details see Ref. J2H|.) The fit only includes data 
satisfying a;, \j\t\ < 160 MeV. Recent work suggests that to describe the backward-angle 
points at 95 MeV requires the inclusion of an explicit delta degree of freedom [531. 

The second error quoted on at N and fix is large because, as alluded to above, cross 
sections depend on the choice of \ijj) at the 10-20% level. This change in the yd cross 
section has a sizable impact on the extraction of and /hv- While not all of the 
wave functions we employed are consistent with yPT, they do all have the correct long¬ 
distance behavior. The differences we see are thus a consequence of different short-range 
behavior in, say, the Nijm93 and NLO yPT NN potentials. As such they should be 
renormalized by 'yNN contact operators. However, Weinberg power counting predicts 
that the operators which do this do not appear until (at least) 0(e 2 P 3 ). Thus the 
degree of variability seen in our results with different wave functions could be viewed as 
inconsistent with this power counting. Indeed, it may be necessary to modify the power 
counting so that, e.g. 7 NN contact operators appear at a lower order than is indicated 
by the naive yPT power counting we have used for (). The issue here is one of relevance 
to nuclear-structure physics: it speaks of the need to properly renormalize the effective 
7 NN operators found when high-energy degrees of freedom are integrated out of our 
EFT of the NN system. Further understanding of this issue is also crucial if \PT is to 
be an accurate calculational tool for low-energy reactions on deuterium. 

5. I thought you said electro weakl 

There is nothing in what I have presented so far that is particular to the photon 
as a probe. Analogous expansions can be made for axial currents. The result is a 
systematic approach to reactions such as pp —> d + e + + u e and ud —> ud. I do not have 
room to properly describe the beautiful work on these reactions by Park, Ando, and 
collaborators BBIBZl. In my view a particular strength of that work is the application 
of the same yPT electroweak operators to A = 2, 3, and 4. Ref. [20] uses the well- 
known tritium beta decay rate to fix the strength of an (otherwise poorly constrained) 
two-body contribution to the Gamow-Teller transition operator. This renormalizes the 
effective NN Gamow-Teller operator to such high accuracy that the pp-fusion S factor 
has an estimated error from higher-order terms of less than 0.5%. The same one- and 
two-nucleon operators are then used to predict the reaction p + 3 He —> 4 He + e + +u e [22] 
and ud scattering E3- This is a great example of yPT’s ability to give a systematic 
derivation of electroweak operators that can be used in a variety of light nuclear systems 
to yield accurate predictions. 
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